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Abstract. 

We review recent extensions of semiclassical multimode laser theory to open 
systems with overlapping resonances and inhomogeneous refractive index. An 
essential ingredient of the theory are biorthogonal quasimodes that describe field 
decay in an open passive system and are used as a basis for lasing modes. We 
discuss applications of the semiclassical theory, as well as other experimental and 
numerical results related to random lasing with mode competition. 
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1. Introduction 

The term "random lasing" encompasses a number of phenomena related to light 
amplification in systems characterized by a spatial distribution of electromagnetic 
field which is much more complex (irregular) than well-defined cavity modes of 
standard lasing structures. Most directly this term describes emission of light by 
spatially inhomogeneous disordered materials not bounded by any artificial mirrors, 
even though it is also used in the case of systems with well-defined cavities, which, 
however, are characterized by chaotic ray dynamics. These two classes of random lasers 
are significantly different, but in some situations, when the statistical properties of the 
field in disordered system and in chaotic cavities are similar, the emission properties 
of both lasers can be discussed on equal footing [1]. In the course of the last decade, 
random lasing has been observed in many different types of disordered materials 
(polymer films [2], porous materials [3], powders [4], ceramics [5,6], clusters [7], 
colloidal solutions of nanoparticles [8]), so that it can be regarded as a universal 
property of disordered structures. 

Lasing in any system is produced due to combination of two factors: optical 
amplification and feedback. In regular lasers the presence of a feedback always assumes 
existence of well-defined phase relations between waves propagating in opposite 
directions. In a different language the presence of a feedback is described as existence 
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of well-defined long-living cavity modes characterized by a regular spatial pattern of 
electromagnetic field, which sets inside the lasing structure in the stationary regime. 
In [9,10], where the concept of random lasers was conceived, a possibility of lasing 
caused by a different type of a feedback was proposed. It was shown that even if 
propagation of light is described in terms of diffusion equation, which completely 
ignores its wave nature and does not have any phase information, one can still have 
a laser-like behavior of emission characterized by a threshold, spectral narrowing, 
relaxation oscillations and other attributes of laser oscillations. Physical origin of 
this phenomenon lies in a significant increase of a length of light trajectory inside 
finite amplifying volume due to multiple scattering. Since amplification results in 
exponential growth of light intensity with the distance traveled inside the gain medium, 
it can be characterized by gain length L g , which depends on the gain factor and the 
diffusion coefficient of light in the medium. The transition to lasing occurs when 
the gain length exceeds the loss length, L;, so that the threshold condition can be 
written down as L g = L\. This type of feedback is called incoherent or nonresonant 
feedback. The latter term refers to the absence of any resonant features in the 
distribution of the field inside the gain medium or, in other words, the absence of any 
distinct modes: the spatial distribution of light intensity is the same regardless of its 
frequency. In this situation the lasing frequency is determined by only one remaining 
resonant element of the system — the atomic transition, so that the single-peak emission 
spectrum is the characteristic feature of lasing with nonresonant feedback [9, 10]. The 
authors of [11,12], where non-resonant feedback was discussed in connection with 
lasing in cavities with rough surfaces, formulated four conditions required to realize 
such a complete nonresonant situations: (i) no mode degeneracy, (ii) equal mode 
loss, (iii) mode overlapping (the resonance width of each mode due to radiative losses 
must be larger than mean intermodc spacing) and (iv) mode mixing (there must be 
processes causing frequency of emitted photons to change by an amount larger than 
the spectral distance between the modes). While conditions (i), (iii) and (iv) are 
realized automatically in disordered systems of large enough dimensions and chaotic 
cavities, the condition (ii) is more difficult to fulfill. Indeed, distribution of widths of 
scattering resonances in chaotic cavities and disordered systems was found to be rather 
broad with different modes having decay rates differing by the orders of magnitude [13]. 
However, when interest to random lasers was renewed in 1994 following the observation 
of lasing from solution of dye molecules surrounded by titanium dioxide particles [14], 
the importance of this condition was apparently underappreciated. Namely, results 
of [14] and subsequent experiments were usually explained on the basis of the concept 
of the nonresonant feedback with its underlying assumption that the transport of 
light can be described within diffusion approximation (see, e.g., [15-19]). Validity 
of diffusion approximation is usually associated with weak scattering of light, when 
the system is far away from Anderson localization transition and interference effects 
are weak. Therefore, nonresonant or incoherent feedback was by extension associated 
with the regime of weak scattering of light (see, e.g., [20,21]). 

Inducing stronger light scattering by increasing concentration of scatterers, it 
became possible to observe a qualitatively new phenomenon [4, 22-24]: with increasing 
scattering multiple emission lines appeared in the spectrum, instead of just a single 
peak at weaker scattering. New peaks, characterized by much narrower line widths, 
emerged one by one with increasing pumping. Since incoherent feedback can only 
result in a single-frequency emission spectrum, it was suggested in [4] that the changes 
in the emission spectrum are due to transition from incoherent to the coherent 
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feedback. This idea was supported by studies of photon statistics, which showed that, 
similar to regular lasers, light emitted at the peak frequencies had Poisson photon 
count distribution [2,25]. 

Originally it was suggested in [4] that the feedback is provided by randomly 
occurring closed trajectories formed by multiply scattered light. Eventually this idea 
was developed to a more general concept of random cavities (resonators), arising 
in strongly scattering medium. Anderson localization was considered as one of the 
possible mechanisms responsible for formation of such cavities [26-28]. In order 
to verify this assumption a great deal of efforts have been devoted to analyzing 
lasing in one-dimensional models in which all states are localized [29-33]. The 
localization-based approach allowed to explain a number of experimentally observed 
results such as mode repulsion and saturation of the number of lasing modes [34]. 
However, it has never been convincingly demonstrated that light in strongly scattering 
three- or two-dimensional samples was indeed close to the Anderson transition. An 
alternative mechanism of formation of cavities that could be responsible for coherent 
feedback in random media was put forward in [35], where it was suggested that 
random fluctuations of the refractive index of a disordered medium can result in 
macroscopically large ring-like configurations capable of trapping light for long times 
and serve, therefore, as random resonators. This model was supported by studies 
of lasing in 7r-conjugated polymers [36-38], where certain degree of universality in 
spectral distribution of lasing modes was found. This universality was explained 
by noting that among multiple random resonators only those with largest Q-factors 
(optimal resonators) mostly contribute to lasing. While distribution of all resonators 
is very broad, the optimal resonators have almost identical characteristics, thus 
explaining observed universality. These random light-trapping configurations are 
analogous to so called prelocalized states, known to exist in the case of electrons 
in random potentials, and can arise even when a disordered system is far away 
from the localization transition. It was, however, shown in [35] that the trapping 
configurations can appear with any appreciable probability only if spatial fluctuations 
of the refractive index are correlated over large enough distances, which might be a 
reasonable assumption for the polymer samples studied in [36-38] , but is more difficult 
to justify for the ZnO powders [4]. A yet another alternative model of random lasing 
was proposed in [39,40], where random lasers were treated as lasers with distributed 
feedback. This approach is somewhat similar to the random-resonator model with 
only difference that, instead of dealing with ring-like resonators, the distributcd- 
feedback model assumed that long-scale almost periodic Bragg-like configurations are 
responsible for lasing. This model, however, has not yet been sufficiently developed. 

More recent developments in the field of random lasers have been associated with 
renewed attention to the weakly scattering samples. Emission spectrum containing 
very sharp spikes was observed in very weakly scattering samples in [41]. It was 
argued in that paper that the appearance of the observed peaks did not require any 
feedback and could be explained by assuming that some of spontaneously emitted 
"photons" travel much-longer-than-average distances. This idea was supported by 
simulations based on random walk model (hence the quotation marks in the "photon"), 
which reproduced fairly well the experimental data. However, it was found in [42,43] 
that there are two types of lines in the emission spectrum of a random laser. One, 
called "spikes" in [42,43], similar to that observed in [41], was detected (with strong 
enough pumping) even in samples without any scatterers at all. The second type, 
called "peaks", only appeared in the presence of scatterers. The spikes and peaks 
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demonstrated significantly different statistical properties, which allowed the authors 
of [42,43] to attribute the origin of the former to amplified spontaneous emission 
and claim that only the latter correspond to true lasing with coherent feedback. 
Significance of these development consists in realization of the fact that weak scattering 
by itself does not guarantee nonresonant feedback and that diffusion model may 
not be applicable to active systems even when scattering is weak. This point was 
reinforced in [27] , where it was shown that even in a weakly scattering active medium 
characterized by strong radiative leakage, as well as spectral and spatial overlap of 
the modes, lasing from individual modes can still occur. Clearly, this is only possible 
because even in a weakly scattering system there exists a broad distribution of radiative 
lifetimes emphasised by the onset of lasing. 

Since observation of multipeak emission spectrum in weakly scattering systems 
often required enhanced pumping and its concentration within a rather small volume 
of the sample, effects of the inhomogeneity of pumping on lasing spectrum have also 
been studied in [8,42]. There it was concluded that the inhomogeneity might play a 
role in promoting this phenomenon. 

As a result of these developments the focus of theoretical research in the field of 
random lasers has shifted from identifying special configurations responsible for lasing 
to accepting that a random laser is a multimode system and needs to be treated within 
the framework of a complete multimode lasing theory. This theory can be developed 
along the lines of standard semiclassical lasing theory, but it has to incorporate such 
features specific for random lasers as a much larger role of radiative leakage of the 
modes and irregular spatial dependence of the refractive index. While this theory is 
far from being complete, a review of recent developments in this area appears to be 
useful. To present such a summary is the main objective of this work, whose structure 
is as follows. In section 2 we discuss several alternative ways to introduce modal 
description of strongly open systems. After presenting the formal multimode theory 
in section 3, we turn to discussion of recent experimental and numerical results in 
section 4. 



2. Modes of open systems 

A problem of introducing modal description for open systems has a long history. 
This problem is important not only for physics of lasers, but also as a first necessary 
step toward quantizing electromagnetic field in open resonators. It is not surprising, 
therefore, that there have been many alternative attempts to introduce the system 
of modes suitable for separation time and coordinate dependence of various physical 
quantities such as electric or magnetic fields. The difficulty of the problem stems 
from the fact that openness makes the problem non-Hcrmitian. Therefore, the 
standard recipes for introduction of modes and quantization based on eigenvectors 
of Hcrmitian operators are not applicable in this situation. In this section we review a 
few alternative methods of defining electromagnetic modes of open systems suggested 
by various authors. 

In the Coulomb gauge V • [e(r) E(r, t)] — 0, electric field E(r, t) is governed by a 
wave equation 

e (r)J^E + Vx (V x E) = 0, (1) 

where Gaussian units with the velocity of light in vacuum c = 1 are used. Often, for 
the sake of simplicity, the vector nature of electromagnetic field is neglected, which is 
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possible if coupling between various polarizations in an open spatially inhomogeneous 
system is insignificant. Then electric field is described by the scalar wave equation 

e(r)^E-^E = 0. (2) 

In this section we consider unloaded, or passive, systems (without gain) and assume 
that all the decay of the field is due to the openness, so that the dielectric constant 
e(r) is real. In order to demonstrate the variety of approaches developed to deal 
with this problem we will discuss in this section Fox-Li modes, subsection 2.1, 
quasimodes, subsection 2.2, system-and-bath approach, subsection 2.3, and the most 
recent development, constant-flux modes, will be considered in subsection 2.4. 



2.1. Transverse Fox- Li modes 

Random lasers are, of course, not the first type of systems, in which radiative 
losses cannot be neglected. Historically earliest, the problem of introducing modal 
description in the presence of strong radiative losses arose in connection with lasing 
properties of so called unstable resonators [44]. Fox and Li [45] suggested to define a 
mode of such a system as a field distribution which reproduces itself after the wave 
makes one complete round trip inside the resonator. One assumes that there is a 
well-defined propagation direction of the wave inside the resonator, so that a solution 
of (2) can be written in the form 

E(r, t) = Re {E(r) exp [i (kz - cot)}} , (3) 

where k — ui\/e [e(r) = const] and E(r) changes on a scale much longer than fc -1 . The 
distribution of the field in the transverse direction, E(r), is characterized by solutions 
of integral equation 



/ 



K(r±,r' ± ,z) il) n (r' ± , z) dr' ± = \ n ip n (r ± , z), (4) 



which formally expresses the idea of the field reproducibility and whose kernel is 
determined by optical characteristics of the resonator. This equation has a form of an 
eigenvector equation for a non-Hcrmitian linear integral operator, whose eigenvalue 
|A„| < 1 describes radiative losses after a round trip. Since the respective eigenvectors 
are not orthogonal, one has to introduce an adjoint operator describing propagation 
in the backward direction as: 

f K(r' ± ,r±, z) <^(r' ± , z) dr' ± = A„0„(r ± , z). (5) 
The two families of functions, called Fox-Li modes, are biorthogonal, i.e., 

t (r±,z) Vn(r_L,^)dr_L = S mn , (6) 
and can be used to construct a modal expansion of the field 

£(r)=^c n («)Wn,z). (7) 

n 

if v n are normalized to unity, then 

K n = J |0„(r ± ,z)| 2 dr ± > 1. (8) 

It was shown [44, 46] that K n is the Petermann excess-noise factor [47] , which describes 
the effect of the openness of the resonator on the fundamental Schawlow-Townes 
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linewidth. While the Fox-Li modes played an important role in understanding unstable 
resonators, their application to random lasers and especially to the problem of field 
quantization is rather limited. 

2. 2. Quasimodes 

The idea of quasimodes was transferred to optics from quantum physics. In order 
to describe scattering resonances in atomic and molecular physics, it was proposed 
to solve Schrodingcr equation with boundary conditions at infinity containing only 
outgoing waves and no incoming incident waves (so called Siegert or Gamow boundary 
conditions) [48] . The solutions of the resulting non-Hermitian eigenvector problem are 
characterized by complex eigenvalues and eigenvectors diverging at infinity. The latter 
circumstance makes it impossible to use these modes as a basis for representation of 
the lasing field or to develop quantization procedure representing field in the entire 
space. At the same time, it was argued in [49] that under certain realistic conditions 
it is possible to use these modes to represent field inside the resonator. For both 
lasing and quantization problem the field is needed everywhere, therefore there have 
been several attempts to reformulate this problem in a way which would produce a 
meaningful system of basis vectors valid in the whole space. 

One of the approaches, developed in [50], uses modes obeying Siegert-Gamow 
boundary conditions to present the field inside the cavity, but uses a different set of 
functions to describe outside field in order to avoid the divergency problem. These 
authors argue that such modes describe "natural" evolution of the field inside the 
cavity after it was created and allowed to evolve freely and, therefore, call them natural 
modes. The authors employ an original definition of the biorthogonal inner product 
by presenting left- and right-propagating components of these modes as components 
of a two-dimensional spinor. This construction was studied in more details for a one- 
dimensional cavity of constant refractive index open at one side and having a perfect 
mirror at the other side. Its eigenfunctions satisfying the outgoing-wave boundary 
conditions are of the form i() n (x) oc e 1KriX + rc~ 1KnX , where n n is complex. Its adjoint 
function (j) n {x) = [ip n (x)]* matches an incoming wave. The corresponding natural 
modes are then 



Natural modes for the external region are defined to have a real wavenumber k, but, 
if continued inside the cavity, they would not satisfy the Dirichlct condition at the 
mirror. Both internal and external modes are biorthogonal to their adjoints in the 
respective regions of space. 

Since the natural modes form complete sets, they can be used to quantize the field 
in the whole space. Namely, the amplitudes of the modes ^ n {x) and <J> n (x) become 
cavity operators a n and b n , respectively. Similar, external operators a(k) and b(k) are 
defined. The system Hamiltonian 



is a sum of internal and external contributions without cross-terms. Coupling between 
internal and external waves in this formalism arises due to noncommutativity of 
the internal and external operators. This circumstance makes application of these 
modes not very convenient, therefore, it would be interesting to try to introduce new 
commuting internal and external operators. This would result in a Hamiltonian where 
coupling would enter explicitly in the standard form of cross-terms. 




(9) 



H = H iD ({a n , b n }) + H out ({a(k),b(k)}) 



(10) 
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2.3. Feshbach projection in the system- and-bath quantization scheme 

Field quantization based on quasimodes [50] (subsection 2.2) is an example of the 
system-and-bath approach: one introduces separate eigenmodes and field operators 
inside the open region designated as the "system" and in the surrounding free space 
("bath"). This procedure can be contrasted with the modes-of-the-universe approach, 
in which the modes are defined in the whole space. The separation into a system and a 
bath often provides a more clear physical description. For example, the complex energy 
of a resonator state immediately yields its lifetime, while extracting this information 
from the real continuous spectrum is less straightforward. 

Feshbach projection technique [51] offers a convenient way to perform the system- 
and-bath quantization in a rather general setting of equation (1) [52], the only 
assumption being that e(r) = 1 outside of a finite domain Q of arbitrary shape. 
The idea of the method is to project the Hilbert space of the modes of the universe 
into the Hilbert spaces of Q and its exterior V with appropriate boundary conditions. 

Separating variables in (1) with an ansatz E(r, t) = Re [E w (r) cxp(— iu>t)] and 
introducing a new, in general, vector- valued function </> w (r) = \A(r) E w (r) we can 
formulate an eigenvalue problem 



Lc(> u = — V x 



= ^cp u (11) 



/e(r)_ 

for the Hermitian operator L. After the projection into the subspaccs Q and V, an 
equivalent system of equations 

Lq Q Lq V \( ^) =co 2( h„\ (12) 

L PQ L-p-p J\v u J \ v u J 

for the restrictions (i^ = 4>u\q an d = 4>J\ V is obtained. The differential operators 
Lqq and L-p-p act in their domains and have the same bulk terms as L (11) and 
interface terms; Lqp and L-pQ act at the interface between Q and V. There is a 
certain freedom in the distribution of the interface terms between the operators, as 
well as in the division of the whole space between Q and V. It is expected, however, 
that physically relevant quantities are independent of this choice [53] . The eigenmodes 
of Lqq and L-p-p are determined from the equations 

i QQ/ i A :=w A/ x A> L vv v m (u) = UJ 2 V m (uj), (13) 

where u\ is real and discrete, u> is real and continuous and m is a discrete channel index. 
The modes are required to obey the boundary conditions that make the interface 
terms of the operators vanish. For example, in one dimension these could be the 
Dirichlet/Neumann conditions on the Q/V side of the interface or vice versa. The 
eigenfunctions /i>(r) and u m (u>;r) of the Hermitian operators form complete sets in 
their domains. Because of the special boundary conditions, an expansion of a mode-of- 
the-universe 4> u in terms of these eigenmodes should deviate strongly from the exact 
result in a layer around the interface. The width of the layer should go to zero as the 
number of terms in the expansions becomes infinite. 

The field is quantized in each domain separately by assigning the annihilation 
operators a\ and b m (u) to the eigenmodes fi x and u m {uj). The system-and-bath 
Hamiltonian 



H = ^hw x a[a\ + ^2 dujhLU bri 
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Wx m {uj)a\b m {uj) + Vx m {uj)a x b m {uj)+YL.c. , (14) 



consists of the internal and external contributions and the interaction part with 
W\ m {u) = {H\\LQv\v m (u)) and V\ m (w) = {n* x \L QV \u m (uj)) . In contrast to the 
Hamiltonian (10), here the system and bath operators commute and the interaction 
terms are explicit. From the Heisenberg equations of motion for a x (t) and b m (u);t) 
one obtains the quantum Langevin equation for a x (uj) in the frequency representation: 



i^2[lj5x\> -Slxx'{u)]ax>(u) + Fx(u) = 0, 



(15) 



where F\(u) is the bath noise operator. The non-Hermitian frequency matrix 

n AV (w) = u x 5 xx , - m (Wt) AV (w) - Axy (w), (16) 
(W+) AA , (w) = ^ W Xm {uj) W x , m (uj), (17) 

m 

contains the imaginary damping term and the real frequency shift, often disregarded. 
In the rotating-wave approximation, which is applicable when the damping is much 
smaller than the typical frequencies, the V\ m (u>) contribution can be neglected. 

Scmiclassically, the electric field in the system (in Q) is E w (r) = 
J2 X a x (uj)fj, x (r)/ \/e(r). Representing the matrix (16) via its eigenvalues and 
biorthogonal left and right eigenvectors as Q(u) = J2k \ r k(v))£lk{v){lk(u)\ we can 
construct quasimodes 



</>fc(^;r) - (/x(r)|r fe (w)), E w (r) = ^ E fe (w; r) 



k 



OiM^f^r)/^) (18) 



where Ofc(w) = (lk(ui)\a(u))) and |/x(r)) [|a(w))] is a column of fJ. x (r) [a x (ui)]. These 
modes depend on the real frequency lo imposed by the exterior. With the help 
of (15) and the bath correlation function at zero temperature (F x (ui) Fy (oj r )) = 
(WW^) xx ,(oj) S(uj — lo 1 ) [54] one obtains the (bath averaged) field correlation in the 
mode k [55]: 

(l k (w)\WW^w)\l k (w)) |</> fe (^;r)| 2 



(E fc ( W ;r)-E£( W ';r)) = 



■ Wfe(w)] 2 + Kftw) c(r) 



-S(u) — lo'), 



(19) 



where Wfc(w) = Ref2/c(w) and Kfe(w) = Imfl^w). Thus, the modes Efc(w;r) indeed 
describe a leaking field with the frequency given by the equation lo = Wfc(w) and the 
decay rate K k (u>), as is expected from quasimodes. 



2.4- Constant-flux states 

So called constant-flux (CF) states were introduced in [56] as an attempt to construct 
a system of basis functions suitable for modal expansion of the field inside an open 
resonator, which would also allow easily calculate field outside of the resonator. While 
these modes are, indeed, very convenient for semiclassical lasing theory, it is not clear 
at the present time if they can be used for quantization of the field. The CF modes 
are very similar to the modes [52] of the previous subsection 2.3: in the interior of the 
cavity region Q they satisfy the eigenvalue equation 



= V x 



V x 



(20) 
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where f2 fc is an eigenfrequency and in the exterior domain V the equation for these 
modes takes the form 



where w is a real external spectral parameter, which does not coincide with the 
eigenfrequency ilk- However, equations (20) and (21) are complemented by outgoing- 
wave boundary conditions at infinity and continuity conditions at the boundary of the 
resonator. One can think of lo as of frequency of some external force exciting the modes 
of the cavity, or alternatively, as a spectral parameter of temporal Fourier transform 
used to convert the problem from time to frequency domain. This spectral parameter 
becomes an integration variable when inverse transform back to time domain is 
carried out. Description of the outside modes by equation (21) not only allows to 
introduce nondiverging energy flux, but also ensures the biorthogonality of the internal 
modes. Indeed, electromagnetic boundary conditions are of mixed type and contain 
dependence on the eigenfrequencies. As a result the orthogonality between adjoint 
modes is destroyed, and transition to (21) is a way to restore it. Similar approach 
was described in the case of elastic waves in [57]. It was shown in [56] that modes 
defined this way form a complete set in the domain Q, which can be complemented 
by a biorthogonal set of adjoint modes. The advantage of these states compared to 
quasimodes obeying Siegert-Gamow boundary conditions is that the former remain 
finite at infinity and describe constant flux of energy coming out from the resonator 
(hence, "constant flux" states). There is no need to introduce special outside modes, 
the exterior field is calculated by simply matching outgoing field with the cavity field. 
This is particularly convenient for calculation of light emission from open resonators. 

The CF wavefunctions and their biorthogonal adjoint functions cf> k (Lo;r) provide 
spectral representation of the interior Green-function operator [56] satisfying outgoing 
boundary conditions at infinity, which can be written in the form 



The indices a, a 1 — x,y,z label the polarization of vector-valued field. The Green 
function satisfying the same boundary conditions can be defined in the system-and- 
bath-approach as [52] Gqq(lo) = [lo 2 — L e ff(u;)] _1 , where the differential operator L e s 
has the form 



Thus, tp k {io) are eigenfunctions of L e g(u>) with the eigenvalues VL 2 k (uj). In fact, it was 
shown in [53] for the one-dimensional case, that the interface terms of L b s(lo) disappear 
precisely when the function it acts upon satisfies the outgoing boundary conditions. 
The poles of the Green function, which correspond to scattering resonances, are 
found by analytical continuation of lo into the complex plane such that the equation 
to 2 = H\{lo) is satisfied. At the same time, in the case of CF states, to is always 
real and, if it is fixed by some external conditions (for instance, it can be frequency 
of incident radiation tuned to be in resonance with the cavity mode), the respective 
resonant eigenfrequency of the CF state obeys a different equation: Re [fifc(w)] = to. 

It is interesting to compare CF modes with those obtained in the system-and-bath 
approach described in subsection 2.3. The first step is to write L c g as a matrix in the 
eigenbasis of Lqq [see (13)]: 



Vx Vxt/^H =Lo 2 xl> k {uj) 



(21) 




(22) 



L e ff(co) = Lqq + Lq-p (lo 2 — L-p-p + ie) 1 L-pq. 



(23) 




(24) 
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where Sl is the real diagonal matrix of u\. The difference between the matrix 
tt 2 (ui) (16) and this matrix reads 

n 2 (u;)- L cS (lj) = [Af7(w),O ] + Afi 2 (w), Afi(w) = Sl(u) - ft , (25) 
where the diagonality of Slo was used. This difference is of the second order in the 
small parameters |AO^v /|^a| and \ui\ — u\i\/\u)\\. Therefore, the CF eigenfrequencies 
and eigenfunctions are expected to become close to ftk(uj) and tp k (io) as the above 
ratios decrease. 

2.5. Open resonators in random-matrix theory 

Explicit construction of quasimodes along the lines of subsections 2.2-2.4 in irregular 
systems may require a fair amount of computations. At the same time, one is 
often interested in statistical description of ensembles of such systems. It is well 
known that certain statistical characteristics of chaotic cavities and diffusive media in 
the short- wavelength limit are universal. This means that these properties depend 
only on the symmetries and, possibly, the boundary conditions, but not on the 
details of the spatial distribution of the refractive index. The universality justifies 
modelling eigenfrequencies and eigenfunctions of these systems using appropriate 
statistical ensembles of random matrices, which can be easier to handle analytically 
or numerically than physical systems. A review of random-matrix theory (RMT) in 
open cavities is given in [58]. 

A convenient starting point is the following, rather general representation of the 
scattering matrix of an open system: 

S(lo) = 1 - 2mW\uo - n^)' 1 W, n eS = n -mWWK (26) 

S(u) is an M x M matrix in the channel space. It transforms a column of incoming 
amplitudes into a column of outgoing amplitudes in some basis of the channel states. 
f2 c fr is called effective Hamiltonian in electronic systems. It is analogous to the 
frequency matrix (16). fl a is an N x N Hermitian frequency matrix of the closed 
resonator; its eigenvalues are the eigenfrequencies of the cavity. The N x M matrix W 
describes coupling between the resonator and the open channels. In the limit usually 
considered in RMT, N — > oo, the frequency dependence of W is neglected. 

RMT prescribes that, in order to model a generic chaotic system with time- 
reversal symmetry, the real symmetric matrix fio should be taken from the Gaussian 
orthogonal ensemble. Without loss of generality, the diagonal (off-diagonal) elements 
of Slo are drawn from a normal distribution with zero mean and the variance of 
2/N (1/-/V). In the limit N — > oo the eigenvalues are distributed according to the 
Wigner semicircle law 

g(u>) = ^l-^, -2<c<2, (27) 

where g(w) is normalized to unity. In the absence of direct coupling between the 
channels, matrix S(uj) should become diagonal after the ensemble average. This 
can be achieved simply by taking fixed (for all members of the ensemble) elements 
W nn = yTra/ 77 > for n < M and W nm = otherwise (the limit M <C N is assumed) . 
Thus, the interaction matrix ttWW^ is diagonal with M nonzero elements j n . The 
coupling strength between the resonator and the continuum is characterized by the 
transmission coefficients 

T n = 1 - \(S nn (u;))\ 2 =2 



In + 7n 

2irg(u)) 



(28) 
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In particular, the coupling is the strongest for 7„ = 1. Since the coupling depends on 
the density g(u>), one has to be careful not to mix the statistics from different regions 
of u> or introduce appropriate scaling in order to obtain universal results. 

Quasimodes of an open chaotic resonator can be modelled via the resonances 
of S(lo), which are the eigenmodes of f2 e ff- Statistical properties of the resonant 
wavefunctions will follow from the surmise that values of an cigenfunction of a closed 
chaotic cavity are uncorrelated Gaussian random variables at points spaced more than 
a wavelength apart [59] . 

Several papers extend the calculations [44,46] of the Petermann factor [47] to 
chaotic systems, where, in particular, there is no separation into longitudinal and 
transverse modes. In the case of chaotic cavity with an opening [60,61], the RMT 
effective Hamiltonian £l c s (26) is modified to account for the amplifying medium in 
the resonator by adding a correction ij a /2 to its eigenvalues, where j a > is the 
amplification rate. As 7 a is increased from zero, the eigenvalues shift upwards in the 
complex plane (the gain is assumed to be the same for all relevant frequencies). The 
first eigenvalue touching the real axis defines the lasing threshold. The Petermann 
factor is shown to be 

K = (l\l)(r\r) (29) 

where \l) (|r)) is the left (right) eigenvector of f2 e ff (for 7 a = 0). A supersymmctric 
calculation [61] reveals that (K) scales as the square root of the number of channels in 
the opening. If the system has a time-reversal symmetry, the matrix Q e g is symmetric, 
i.e., W is real. Then one can impose the condition \l) = \r)* , which leads to an 
alternative expression 

^TRS = (r\r) 2 . (30) 

In a chaotic dielectric resonator (a domain of uniform refractive index n > 1 
surrounded by a medium with n = 1) light can leak anywhere along the boundary. 
Hence, the number of output channels M scales with the size N of Q e g, breaking 
the standard RMT assumption M -C N. A combination of RMT with the Fresnel 
laws [62] produces a scattering matrix 

S(u>) = -R + TF(uj)[1- RFiuj^T, (31) 

where F{lo) is the intracavity propagator and R and T are diagonal matrices of 
reflection and transmission. The Petermann factor is of the form (29), but the 
eigenvectors refer to the matrix RF(u). After a symmetrization procedure, the 
form (30) can also be obtained if time-reversal symmetry is present. Comparison 
of the RMT calculations with the quantum-kicked-rotator model points at loss of 
universality as n becomes close to unity, i.e., in a strongly open resonator. 

It should be noted that the near-threshold treatment [60-62] yields the value 
of Petermann factor which is twice as large as the correct result. This happens 
because phase diffusion (which mainly contributes to the line width) is amended at 
the threshold by the amplitude fluctuations. The correct prefactor was inserted "by 
hand" in equations (29) and (30). A linearization of the quantum Langcvin equations 
far above the threshold [63] yields again the result (29), this time with the proper 
prefactor. 
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3. Multimode laser theories for open and irregular systems 

As it was mentioned in introductory section 1, the current trend in theoretical 
description of random lasing consists in extending standard semiclassical mutlimode 
lasing theory [64, 65] to situations specific for random lasing: modes with broad 
distribution of radiative lifetimes and irregular spatial patterns, inhomogeneity of the 
background refractive index, large number of lasing modes emitting in the regime when 
pumping significantly exceeds the threshold value. While it is true that the multimode 
lasing was considered one of the signatures of coherent feedback, the discussion of this 
issue in earlier works related to strongly scattering systems was mostly concerned 
with a rather trivial situation, when multiple modes originated from different non- 
overlapping cavities (localized states), which were assumed to be in the single- mode 
regime due to mode competition [30,32,34]. The only nontrivial multimode effect 
reported in strongly scattering systems was observation of mode coupling in [32], 
however, an attempt of its explanation using a two-mode lasing theory given in [66] 
relied on incorrect form of rate equations (this point will be explained later in this 
section). 

One of the first attempts to apply multimode lasing theory to chaotic lasers 
was undertaken in [1], where equations of standard semiclassical laser theory were 
combined with ideas of random-matrix theory concerning statistical properties of 
eigenmodes of chaotic resonators. Cavities studied in [1] were characterized by a 
relatively long radiative lifetimes of the modes. Statistical properties of random 
lasers with strong radiative losses were studied in [67-69] using combination of 
random-matrix theory with the method of Feshbach projection. A statistical study 
of one-dimensional strongly open random lasers based on calculations of lasing modes 
"from the first principles" rather than on phenomenological random-matrix-type 
considerations was carried out in [70]. 

The role of nonuniformity of the refractive index in formation of lasing modes 
and its consequences for lasing dynamics was discussed in [71,72], where the idea 
that lasing modes can be significantly different from modes of cold cavities and must 
be determined selfconsistently was formulated. Application of the theory [71,72] 
to random lasers was, however, limited because of restriction of the selfconsistency 
requirements to the linear regime only and neglect of radiative losses of the modes. 
Both these limitations were removed in [56, 73, 74] (see also recent review article [75]), 
where spatial structure of lasing modes and lasing frequencies were determined 
selfconsistently from fully nonlinear theory. The approach developed in these papers 
is based on three main ingredients: (i) use of CF states to incorporate radiative losses 
of the system, (ii) neglecting population pulsation, which allowed to take into account 
nonlinear interactions up to infinite order in field intensity and (iii) determination of 
lasing modes and their frequencies selfconsistently. Results of those works revealed the 
presence of strong effects related to nonlinear interaction between modes in systems 
with strongly overlapping (both spectrally and spatially) modes. However, since this 
method is based on numerical computation of fixed points of a certain nonlinear 
map, which, as it often happens in nonlinear systems, might have multiple stable 
and metastable points, the usefulness of this approach might be limited to systems 
with relatively large intermode spacing. Its applicability to truly diffusive random 
lasers is difficult to assess from the data published in [56,73-75], since they do not 
contain any information about mean free path, Z m f p , of light in studied structures. 
Diffusive regime arises only when relations R ^> Z m f p ^> A, where A is the wavelength 
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of light and R is the size of the sample, are satisfied, which requires much larger 
samples than the condition R 3> A, fulfilled in samples studied in [56, 73-75]. 

In this and the following section we derive a semiclassical multimode lasing theory 
taking into account recent achievements discussed above and present some examples 
of its application to situations relevant for random lasers. 



3.1. Semiclassical laser equations 

A starting point for the semiclassical description of lasers [64, 65] is the wave equation 
for the electric field (1) with the polarization P(r, t) as a source that generates the 
field (c = l): 

e (r)^E + Vx(VxE) = -4^P(r ;i ). (32) 

In the simplest model, the polarization is produced by two-level active atoms and 
obeys the equation 

dt 2+2l± dt + l/ J p = ~ 2 ^T E(r ' t] An(r ' t] ' (33) 

where An(r, t) is the population-inversion density, d is the magnitude of the atomic 
dipole matrix element, v is the atomic transition frequency (homogeneous broadening 
is assumed) and jj_ is the polarization decay rate. The population inversion, in turn, 
depends on the electric field and the polarization, 

j t An - 7|| [Ano(r, t) - An] = J-E(r, t) ■ ^P(r, t). (34) 

If the right-hand side vanishes, An relaxes with the rate 7|| to the unsaturated 
population invertion Ano(r, t), which is determined by the pump. 

The coupled equations (32)-(34) yield, in principle, distribution of the electric 
field in the system, if Ano(r, t) is given. Unlike the traditional derivation of lasing 
equations, which is usually done in time domain, we find it more convenient to proceed 
using frequency representation. We introduce the Fourier transforms 

1 r°° 

E(r,t) =-Rc/ dwE w (r)e" llJi , (35) 

JO 

P(r,t) =-Re/ d^P u (r)e- lwt , (36) 
Jo 

An(r,t) = — / do;An u (r)e- lut . (37) 

2n J-oo 

in the form chosen to facilitate application of the rotating-wave approximation. 
In addition, we assume that the time dependence of the field and polarization is 
determined by fast oscillations with frequencies which are close to atomic frequency v 
and residual slow time dependence. In the frequency domain this means that only 
Fourier components E w and P w with u) in the small vicinity of v contribute significantly 
to the dynamics. Therefore, after performing Fourier transform of (32)-(34), one can 
neglect terms of the order of (oj — v) 2 in uo 2 = (w — v + v) 2 w v 2 + 2u{uo — v) and write 
down the resulting equations as 

-e(r) {-v 2 + 2vu) E„|Vx(VxE„) = 47r^ 2 P^, (38) 

d 2 r°° 

[-i(w-!/)+ 7 i]P u = -i— / dw'E^An w _ w ,, (39) 
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(-iw + 7||)An w - 27r 7 ||An (r) ( 5( W ) - — jf <L/(E;,_ W • P w , - E u , +U • P£,), (40) 

where we assumed that the pump Ano(r) is time independent. 

In the linear approximation [71,72] (valid below and not far above the lasing 
threshold) we neglect the quadratic terms in (40) and use 

An£°> = 27rAn (r)%) (41) 

in (39). Polarization in this approximation is given by 

Pl 1} = -i(d 2 /h 1± )D(Lu)An a (r) E w (42) 

and, when substituted to the right-hand side of (38), yields the following equation for 
the electric field: 

d 2 2 

-e(r) {-v 2 + 2vuj) E„|Vx(VxEj = -4iri-^-D(uj) An (r) E w , (43) 

where D(u>) = [1 — i(oj — v)/j±]~ 1 . 

One can use any of the system of modes T/>fc(w; r) discussed in subsections 2.2-2.4 
and their adjoint modes 4> k (u;r) with respective eigenfrequencies (we use the 

same notation for all types of modes) in order to generate modal expansion (18) 
of electric field in equation (43). The expansion coefficients are found using the 
biorthogonal functions as 

a fc (w) - J dr^/I^) < j>l(uJ■, r) • E„(r). (44) 

Then (43) is reduced to the matrix eigenvalue problem 

^2 [uo~kk> - Hfefe/ (u)] a k , (uj) = (45) 



k' 



with 



- d 2 

f2 fc fe'(w) = Vt k {u))5 k k> +\2nv- — D (uj)Vkk' (oj) , (46) 

V kV (w)= J dr0£(u,;r)^ fe ,(u,;r)^M (47) 

The matrix V kk ' (oj) becomes diagonal for uniform Ana and e. The lasing thresholds 
and frequencies at the thresholds are determined from the system of equations 
uj = Keflk(oj) and Imflk(uj) = 0, where Ctk(u>) are eigenvalues of Clkk'{w)- 

3.2. Third-order theory 

Nonlinear effects can be included by iterating equations (39) and (40), with the 
field as a small parameter [64,65]. Namely, Pi, is inserted in (40) to obtain the 
correction Anl, , which, in turn, is used in (39) to yield the contribution to the 
polarization p£ 3 ^ , of the third order in the electric field. 

Assuming that lasing modes exist, we can present field as a sum of oscillating 
terms with slowly varying amplitudes, 



E(r,i) = Re 



^Ej(r,i) e" iW! * 



E ( „(r)=5^E I (w-a; I ;r) ) (48) 
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and positive frequencies u>i close to v, which are to be determined selfconsistently. 
If the amplitudes E;(r,i) vary slowly on the scale of u>i, their Fourier transforms 
E;(w — u>i; r) are strongly peaked at u>i. The slow- varying- amplitude approximation in 
the frequency domain amounts to the replacement E;(w — uji; r) — > 7rE;(r, i) 5(w — wj) 
in nonlinear terms before performing frequency integrals when transforming to the 
time representation. Further, we average out interference terms, oscillating at the beat 
frequencies, and neglect mode degeneracies (|w;—(jj m | 3> |E;|/|E;|). After the nonlinear 
correction p£ 3 ^ is added to the right-hand side of (38) and the linear contribution is 
diagonalized according to (45), we arrive at the third-order lasing equations 



|^ + i \pk{u m ) - w m ] | a km {t) = 



-— (J^J DM |dr^^0;( Wm ;r).£{2Rc[^)]E„ l (r,i)|E i (r,t)| s 



+ (1 - S ml ) D\\(u m - ui) [D{io m ) + D*(un)] E,(r, t) [E*(r, t) • E m (r, t)]} . (49) 

akm is the amplitude of the fc-th component of the m-th lasing mode in the basis of 
eigenfunctions of linearized problem (43) and (45): 

N b 

E m (r,i) = ^2a km (t)tp k {u) m ;r)/^/e{r), m = 1, . . . ,N m . (50) 
fe=i 

This is a system of A^b x N m equations, where is the basis size and N m is the 
number of lasing modes. The terms proportional to D\\ (cu) = (1 — ic<j /t|| ) — 1 arise from 
the population pulsations at the beat frequency cu m — uji interfering in (39) with the 
oscillations at the frequency u>i (and, thus, surviving the averaging). In the stationary 
regime (dafc m /dt = 0), the amplitudes Hkm and frequencies ui m can be determined 
from equations (49) by an iteration procedure (see section 3.3). 

The number of equations in system (49) can be reduced to Ab if one assumes 
that the wavefunctions of the lasing modes are given by the linear approximation (43) 
and do not have to be determined selfconsistently from the nonlinear equations (49). 
This is the approximation used in [71,72], where it was justified by the fact that 
such a partially (in the linear approximation) selfconsistent treatment still allows one 
to eliminate fast oscillating terms in the nonlinear polarization and to introduce the 
slow-varying-amplitudc approximation resulting in rate equations for the respective 
amplitudes. In this linearly selfconsistent approximation the amplitudes can be 
presented as a km (t) = a m (t) 5k m and the lasing equations for a scalar field take the 
form (cf. [67, 69, 70]) 

{d7 +i [^ m ( Wm ) ~ ja m = ^T/^jj" (^[) D ( LJ m)a m ^2B m i \ai\ 2 

x {2Re [D(ui)} + {l-5 m i)D\\(u m -ui)[D(ij m )+D*(iJi)]}, (51) 

B m i=V J dr^^p-^(w m ;r) ^ m {uj m ;r) |^(wj;r)| 2 , (52) 

where V is the volume. Equations (51) generalize the standard third-order 
semiclassical theory with the saturation (hole-burning) terms [64,65] to the case of 
strongly open and irregular systems. Separating real and imaginary parts of this 
equation one obtains rate equations for intensities I m = \a m \ 2 of the modes and an 
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equation for lasing frequencies: 



2ttv ( d 2 
Vhry\\ \hj± 



) 



2 



— - 2 Im [n m (uj m )] [ I m = 



7 m Re D(w m ) ^ (■ ■ ■) 



(53) 



) 



2 



Re [f2 m (w 



'm — 



F/l7|| V^7-L 



Im £>K,) ^ (■ ■ ■) 



(54) 



where the sum ^ z (• • •) appearing in (51) depends on intensities of all lasing modes. 
These rate equations do not contain any linear coupling terms contrary to the 
assumption made in [66]. 

These equations show that all statistical characteristics of laser emission 
(frequency, threshold and intensity distributions) are determined by certain integrals 
involving cigenfunctions of cold cavities. Transition from strong to weak scattering 
manifests itself in changing statistical characteristics of respective quantities. However, 
in spite of large amount of work on wavefunction statistics in closed systems, the 
statistical properties of self- and cross-saturation coefficients in open resonators have 
not yet been studied. At the same time, it is clear now that this statistics is responsible 
for various regimes of behavior of random lasers [70]. We will discuss this point in 
more details in subsection 4.1. 

3.3. All-order nonlinear theory in the time-independent population approximation 

It is possible to obtain lasing equations valid in all orders in the electric field in 
a closed form if one neglects time dependence of the population inversion. As seen 
from (49), the population-pulsation contribution can be neglected if \D\\ (u) m — u>i) | <C 1. 
Typically, the lasing modes are excited within the gain bandwidth 7j_ around the 
atomic frequency. Then, the above condition reduces to 7j_ 3> 7|| . 

Requiring that An w (r) = An(r) S(u>), we express P w from (39) and insert it 
in (40). If this assumption were actually consistent with the equation (40) one 
would, after carrying out mode-of-the-field expansion (48), end up (in the frequency 
representation) with terms, which were also proportional to 5(u>). In reality, in 
addition to "correct" terms one would obtain a number of terms proportional to 
5 functions of various combinations of lasing frequencies, which describe oscillations 
of the population. Neglecting this "oscillatory" terms is equivalent to keeping only 
diagonal contributions |E;| 2 in the modal expansion of the E* • P term, quadratic in 
the field. In this approximation An(r) can be determined selfconsistently and inserted 
into (39) and (38) to obtain the lasing equations 



In contrast to (49), the field here is expanded in the quasimodes of the passive system 
with the frequencies f2fc(cj m ), while the linear mode coupling is included in the right- 
hand side. These equations represent generalization of time-independent equations 
derived in [56], which are obtained from (55) by assuming time independence of the 
respective amplitudes. Equivalcntly, equation (55) can be derived by treating the 
polarization term in (38) as a source and using the Green function to write down the 




(55) 
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solution of this equation as [56]: 

E w (r) = -^ [ dr' C - 1 /2( r ')G(a;;r,r')P w (r') 
\A(r) V 

,_ 7|| E,Re[^)]|Ei(i-,*)l 2 ' (56) 

Replacing the Green function with its spectral representation (22) in the rotating- 
wave approximation and integrating as in (44), we again arrive at equation (55). In 
the stationary case this equation is reduced to a nonlinear eigenmode problem 

^ Tkk'i^m) a-k'm = P^ 1 dkm, (57) 

fryxw-fikM./ e(r') 1 + ^ R e [£>( Wl )] |E,(r)P k ' 

where the unsaturated population inversion Ano(r) = p5no(r) is split into the overall 
pump strength p and the pump profile <5no(r). The field distribution in mode m is 
E TO (r) = J^k a km' l Pk( U} rn] r )/\A( r )- If the basis of constant-flux modes (section 2.4) 
is used, the field outside of the system can be obtained by continuation. 

An algorithm to determine the lasing-modc frequencies oj m and expansion 
coefficients ctkm as p is increased continuously from zero is described in [74]. Below 
the threshold, where all ak m = 0, one looks for the eigenvalues of the linear T(u>). 
Changing u, the eigenvalues can be made real, one at a time. The largest real 
eigenvalue p^ 1 yields the threshold pump strength and the corresponding lu is the 
lasing frequency at the threshold. Above the threshold, for p > pi, the pump is 
increased in small steps and the solution aki{p) for the first mode is determined 
iteratively from (57). The second mode appears when the second largest eigenvalue 
P2 1 of T(uj) linearized "around" the first mode becomes equal to p~ x . The procedure 
is continued to find higher modes. 

4. Examples and properties of multimode random lasers 

4-1. Threshold and number- of -modes statistics 

Distribution of thresholds and average number of lasing modes as a function of pump 
strength was calculated in [1] for an ensemble of weakly open chaotic cavities. Each 
cavity was opened via M small holes (diameter <§C wavelength) , which together carry 
M open channels. The passsive-mode decay rates k in this system are distributed 
according to the x 2 distribution with M degrees of freedom, 

p M= {J ^^y M,2 - 1 ^[-^y), v = k/W, (59) 

where T(x) is the gamma function. The distribution is wide for small M and, for 
M = 1, it increases as y~ x l 2 when y — > 0. This property leads to a wide distribution of 
lasing thresholds, Ano,thr, which behaves as An^f^ 1 for small Ano.thr- The average 
threshold is much less than the nominal value An — (k) hj±e/2Trvd 2 , which is the 
pumping required to overcome the average loss at u> — v [cf. (46)]. An increase of 
threshold fluctuations with localization was observed in a one-dimensional disordered 
model [42]. 
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Considering the rate equations (53) in the stationary regime and neglecting the 
population-pulsation term containing Dn(u> m —u>i), one obtains a matrix equation for 
the intensities I m — \a m \ 2 of lasing modes: 



^ An y m 2d 2 V C t 
> A mX Ii=\-- — — , A„ d = — —Re 

An C m e/i 2 7_L7|| L m 



D(w m ) B n 



V 2 An 



(60) 



where y m = n m /{n) (n m is the decay rate of mode m) and C m = ReD(w m ). This 
equation must be complemented by the condition that it has only positive solutions. 
The rest of the basis are nonlasing modes and their intensities are set to zero. It 
should be noted, however, that the positiveness of lasing intensity does not, by itself, 
guarantee that the found solution is stable, the fact well known for simple two- 
mode models [64]. The stability of the solutions can only be verified from the time- 
dependent equations (53), therefore, the estimates of the number of modes based on 
time-independent equations cannot, in general, be considered as completely accurate. 

For a weakly open chaotic cavity [1] one can assume that (a) the cigenfunctions 
are almost real and orthogonal and (b) they can be described as random Gaussian 
functions which are uncorrelated for different modes [59]. If, in addition, one assumes 
that the background dielectric constant, e(r), and pumping rate, Ano(r), are both 
uniform, the correlator (52) takes the form of 

B m ie 2 /V 2 An = 1 + 2<J mi . (61) 

In this case matrix A m i can be inverted analytically yielding a dependence of mode 
intensities on the pump strength. Let us order the modes m = 1,2,... in the order 
they are excited as the pump increases. It can be shown that, in this case, ymj^m 
form an increasing sequence. A threshold condition for the mode m results in the 
equation [1]: 

fm ,,\Vm if yi An 

l2 +1 Jz^-2gzr av (62) 

which relates the number of excited modes iV m = m and the pump strength An in 
a particular cavity. The ensemble average (iV m ) is calculated [1] with the help of the 
distribution Pniy) for y <C 1 and is found to scale as Ari^' M+2 ' asymptotically for 
large pump strength. 

In [67] (iV m ) was computed for nonweakly open cavities modelled by random 
matrices. In particular, validity of the relation (61) was studied numerically for strong 
coupling to the bath -f n = 1 (section 2.5). It was shown that for n < (k) the assumption 
of uncorrelated eigenvectors works well, while for larger k deviations from Gaussian 
statistics become stronger. 

The power-law asymptotics for (N m ) was confirmed by numerical simulation of 
decay rates entering (62) using random matrices with )„ « 1 [68]. The standard 
deviation crjv m varied as An^^ M+2 \ In the case 7„ = 1, the ratio <rjv m /(A^ m ), but 
not <7iv m and (N m ) separately, obeyed a power law with an exponent that depended 
on M. The difference in the results for the weak and strong coupling occurs, because 
the decay-rate distribution for j n ~ 1 is no longer of x 2 type [58] . 

In a one-dimensional disordered system [70] assumption (61) fails completely, 
indicating a different type of mode competition compared to chaotic systems. 
Numerical simulations showed that the number of lasing modes saturates below the 
basis size with increasing pump. This effect is related to the nonmonotonic dependence 
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of mode intensities on the pump strength and complete disappearance of some modes 
for pump exceeding certain thresholds. 

The mode suppression was also reported in disordered disk lasers studied within 
the theory of section 3.3 [74]. Comparing the dependences of lasing frequencies and 
intensities on the pump, it was noticed that when two modes come close together in 
frequency, one of them can be suppressed. The mode thresholds and intensities, but 
not the frequencies, were found to be very sensitive to the pump spatial profile. 

Jf.,2. Frequency and intensity statistics 

It is well known that passive closed chaotic systems without spatial symmetries have 
level repulsion, i.e., the probability density for zero frequency spacing vanishes. A 
natural question is, how the spacing distribution for lasing modes in a random laser 
is connected to the distribution for passive modes in the underlying system without 
gain? It can be seen in the following examples that mode selection and competition 
normally enhances the repulsion in a laser. When lasing modes are close to passive 
modes, this property is rather obvious: even if two passive modes cross, not both of 
them will necessarily lase. 

Spacing distributions in two-mode chaotic lasers modelled with random matrices 
were computed numerically in [69]. Mode repulsion was present both in the cases 
of weak and intermediate openness, 7„ C 1 and 7„ = 1 (section 2.5), even though 
the passive frequencies can cross for 7„ = 1. (If two passive modes have the same 
frequencies, they have quite different lifetimes due to repulsion in the complex plane. 
Hence, only one of the two modes will be lasing for moderate pump strength.) When 
the gain bandwidth j± is close to the mean level spacing Aw in the passive system, the 
spacing distribution for the lasing modes is well described by the Wigncr surmise [76] 



derived for passive closed chaotic cavities in the random-matrix theory. Again, 
this form works also for 7„ = 1, when the spacing statistics for passive modes is 
closer to Poissonian. This example shows that formal coincidence between a spacing 
distribution for lasing modes and the Wigner surmise does guarantee that the physics 
behind them is the same. 

Mode repulsion with deviations from the Wigner surmise was found numerically in 
one-dimensional disordered lasers [70] . There were two reasons for the repulsion. First, 
some modes were coupled, because the system was not very long and the modes could 
overlap. Second, when two modes were localized (and could have close frequencies), 
the mode that was closer to the opening had higher threshold and was not excited. 

When two modes have close frequencies, it may become necessary to take into 
account the dependence of mode frequencies on the pump. As mentioned above, a 
correlation between mode repulsion and suppression during the change of the pump 
level was observed in a numerical study of two-dimensional disordered laser [74] . 

Spacing distributions were measured in colloidal solutions containing TiC>2 
scatterers and a laser dye [43]. The system was in the weak-scattering regime 
in the sense that the scattering mean free path was much longer than the pump 
excitation cone. The lasing frequencies were more or less regularly spaced, exhibiting 
the mode repulsion. The spacing distribution had a maximum, but could not be 
fitted well with the Wigncr surmise. The average mode spacing scaled with the 
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dye concentration. (Increasing the concentration reduced the gain volume, which 
led to the reduction of the number of modes.) The spacing fluctuations increased 
with the scatterer concentration. Some of the experimental results were supported by 
numerical simulations in a onc-dimensional disordered system at the threshold. The 
statistics of lasing peaks was compared with the statistics of spontaneous-emission 
spikes that appear in the background of the emission spectra. The spikes were 
attributed to photons created in single spontaneous-emission events and amplified 
over long paths [41,77]. Coherent feedback is not required for the appearance of 
spikes. The spikes' positions in the spectrum were uncorrelated, which was reflected 
in the Poisson spacing statistics. 

Spectra with almost equally spaced lasing peaks were obtained in Ti02 colloidal 
solutions with strong reabsorption outside of the pumped volume [8]. Weak scattering, 
on the one hand, and reabsorption, on the other hand, result in an effective cavity 
being formed by just two scatterers located at the ends of the excitation cone. (The 
"cavity" has the maximal possible length, because the gain grows exponentially with 
the path length, while the probability of photons leaving the cavity scales as a power 
of the length.) The effective cavity is of the Fabry-Perot type, therefore the lasing 
peaks are equidistant and the spacing scales inversely with the cone length. 

Another system that shows mode repulsion is porous GaP filled with dye 
solution [3]. The transport mean free path was of the order of A and about three 
times smaller than the pump spot size. The spacing distribution could be roughly 
fitted with the Wigner surmise. 

The difference between lasing peaks and spontaneous-emission spikes in TiC>2 
colloidal solutions with weak reabsorption [43] (see above) emerges also in the emission 
intensity statistics. Two statistical ensembles were considered: (1) intensities collected 
from all wavelengths in some range /(A)/ (/(A)), normalized by the intensity averaged 
over many shots at given A, and (2) peak and spike heights of the functions /(A) / (/(A)) . 
For both ensembles the probability distribution showed similar asymptotic behaviour 
at large intensities: it had a power-law tail above the lasing threshold and decayed 
exponentially below the threshold or in the absence of scatterers (neat dye solution). 
Numerically computed distribution of lasing-modc intensities in a one-dimensional 
disordered laser [70] had a power-law decay, as well. 

It should be mentioned that a power-law asymptotics may also appear in lasers 
with incoherent feedback near the threshold [78]. Thus, caution should be exercised 
when using intensity distribution as a test for coherent lasing. 

4-3. Structure of lasing modes 

One of the important recent development, which is relevant not only for random lasers, 
but for the entire field of laser physics is the realization of the fact that so called lasing 
modes may differ significantly from modes of passive cavities. It was noted in [71] that 
spatial nonuniformity of the refractive index and pumping can result in gain-induced 
linear coupling between modes of passive cavities, which results in formation of new 
modes. These ideas were taken further in self-consistent theory of [56, 73-75], where no 
a priori assumptions about the spatial structure of lasing modes were made and they 
were found from nonlinear self-consistent equations (55). Calculations of [56,73-75] 
found significant modifications of the spatial profile of lasing modes due to nonlinear 
mode coupling. However, as we already mentioned, it is not clear if the systems studied 
in [56, 73-75] can be considered as being in diffusive regime. At the same time, no 
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changes in the spatial structure of a lasing mode with increased pumping intensity 
were found in [27,79], where the studied structures were clearly identified as being 
either in localized [27] or diffusive [79] regimes. 

Modification of modes due to presence of gain was observed in numerical 
simulations of a one-dimensional random laser below or at the threshold [42,43], 
but only in the presence of spatially nonuniform (local) pumping. While modes in 
a uniformly pumped system were close to passive modes and their intensity grew 
exponentially towards the boundary, the transition to local pumping changed them 
substantially: they did not grow exponentially outside of the gain volume, but still 
were extended over the whole system. The number of lasing modes under local 
excitation was found to be less than the number of passive modes in the same frequency 
range, but larger than the number of passive modes in the reduced system defined by 
the pump region. In the case of nonuniform pumping the mode modification appears 
already in the linear approximation [72], so that the roles of nonlinear effects in this 
simulation is unclear. 

Structure of lasing modes was also studied in a system with weak scattering and 
strong reabsorption [8] . It was found that in such systems the gain volume surrounded 
by strongly absorbing medium forms an effective cavity, where lasing modes are 
localized. Numerical simulations below and at the threshold in two dimensions showed 
that the lasing modes in this case are very close to the passive modes of the effective 
cavity. 

5. Conclusions 

Current research on multimode random lasing is moving along several major directions, 
such as (i) extension of conventional laser theories to open and irregular systems; 
(ii) statistical properties of lasing modes; (iii) mechanisms of random lasing (quantum 
and classical localization, extended modes), to name a few. Up to date a large number 
of experimental and numerical observations are collected. However, the systems 
are hard to access analytically, as they consist of a number of strongly interacting 
components (electromagnetic field, gain medium, scatterers, boundaries) and several 
factors (openness, disorder, nonlinearity, noise) are not negligible within a wide range 
of parameters. 

Recent important developments in the semiclassical multimode theory and 
random-matrix theory added to understanding of the properties of lasing modes in the 
stationary regime. To enable a direct comparison with experiments, mostly performed 
under pulsed-pumping conditions, it would be desirable to study time-dependent 
behaviour and relaxation processes. A detailed analysis of mode stability, hysteresis 
phenomena and quantum-noise effects in the stationary regime is also lacking. To 
address the role of localization in lasing feedback at an adequate level, numerical 
simulations of more realistic (three-dimensional) models might be necessary. 
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